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In the standard procedure for calculating the decay rate of a metastable vacuum the solution of the 
classical Euclidean equation of motion of the background field is needed. On the other hand radiative 
corrections have to be taken into account already in the equation of motion. Hence, the latter one 
has to be the functional derivative of the effective action with respect to the background field. This 
is of crucial importance in theories in which the symmetry breaking is due to radiative corrections. 
Usually the effective potential is considered only, neglecting the corrections due to the derivative 
terms of the effective action. In this article a bounce solution from an equation of motion which takes 
into account the full effective action in the one-loop approximation is calculated. A computational 
method that yields a strict separation of the divergent contributions to the effective action from 
the convergent ones is obtained. This allows a wide freedom in the choice of regularization and 
renormalization schemes. The model under consideration is the SU (2)-Higgs model. The fluctuations 
of the complete bosonic sector, i. e. gauge field, Higgs and Goldstone boson contributions, are 
taken into account. The bounce is then self-consistent to one-loop order. The obtained results for 
characteristic quantities of the transition as the nucleation rate and the number of nucleated bubbles 
per volume are compared to other, non-self-consistent approaches. 

ll.15.Kc, ll.lO.Gh, ll.lO.Wx, 98.80.Cq 


I. INTRODUCTION 

When dealing with field theory in cosmological context, one often encounters the problem that a field is not in a 
state which is the absolute minimum of the potential, the true vacuum, but in a metastable local minimum, the false 
vacuum, which is at higher free energy. The transition from false to true vacuum proceeds by a quantum mechanical 
tunnelling process or by a classical transition induced by thermal fluctuations for a system at finite temperature. The 
transition is local in space. During the transition a bubble is formed which is a region in space with true vacuum 
surrounded by false vacuum. The bubble just large enough not to collapse is called critical bubble. Once nucleated, 
the bubble expands, converting false vacuum to true vacuum. 

If the mass of the Higgs boson is not too large the electroweak phase transition is of first order and proceeds by 
bubble nucleation. This scenario is interesting because of the possibility to explain the baryon asymmetry within the 
minimal standard model [Qj^. 

A quantity of fundamental importance for a phase transition is the transition rate. The nucleation rate per volume, 
7 , can be calculated from the bounce-solution of the Euclidean field equation. This solution is also called classical 
solution. For a theory with only one scalar field $ with a potential U ($) one has to solve the equation of motion: 

- ^C/(d>) = 0 . (1) 

The boundary condition is, that the solution ^ = (j) tends to its false vacuum value if one of the tends to cx). The 
nucleation rate can be written in the form 


7 = Aexp(—H) , 

where B is the Euclidean action of the classical solution (fi: 


( 2 ) 
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B = 


( 3 ) 


+ U[4>) 


In the semiclassical approximation the factor A is given by an expression involving functional determinants [|5|-|7[] . The 
negative logarithm is the one-loop contribution to the effective action. 

If U{^) only has a single minimum the bounce-solution does not exist. Nevertheless, in some models the vacuum 
structure is not determined by U (<&) alone but it is changed by radiative corrections. These can be calculated by taking 
into account the quantum fluctuations of the scalar field and those fields coupled to the scalar one §. Depending 
on the specific model, the radiative corrections yield an effective potential with different minima [p|-p^. The 

bounce-solution is then calculated from a held equation that contains this effective potential. However, the effective 
potential is only a part of the effective action. It does not take into account the space-time dependence of the 
background fields, i. e. terms, that depend on derivatives of the background fields. 

This raises the question, whether it is justified to neglect these terms in calculating the bounce ||^. In case of 
the electroweak phase transition the bosonic as well as the fermionic contributions to the effective action have been 
calculated, based on a bounce-solution, that was determined from the effective potential and not the effective action 
00- In im an improvement to take into account non-local effects already for the calculation of the bounce was done 
by modifying the kinetic term of the background field. The radiative corrections were calculated using approximations 
as the derivative expansion |l^ and the heat kernel method Jl4[ or by exact numerical computation schemes 00- 
The corrections due to the bosonic sector (gauge and Higgs fields) were found to be of the same order as the classical 
action. In such a case the non-local terms of the effective action are not negligible and should be taken into account 
when calculating the classical solution. As the fermionic corrections are small compared to the classical action [p^ , 
their non-local contributions to the effective action can be omitted in the equation of motion. 

The standard procedure to calculate the nucleation rate works in the following way: First one adds the one-loop 
corrections of the effective potential to the classical action B and gets an approximation S of the one-loop effective 
action. Based upon the functional form of S the field equation is derived and the bounce is determined. This is used 
to evaluate S and the prefactor A of the transition rate which is given in more detail by 


7 = 


LO- 

27r 



exp ( -S' - - In Jb 


( 4 ) 


Here is the modulus of the eigenvalue of the bubble’s unstable mode. Jb denotes the product of all fluctuation 
determinants of the bosonic sector. Aside the dominant fermionic contribution due to the top-quark which is usually 
taken into account for the setup of S the fluctuation determinant due to fermions is neglected here. The aim of this 
paper is to improve the calculation of the bounce which is then used as the numerical input for the evaluation of the 
nucleation rate. It is no longer determined as a saddle point of S but as a saddle point of the full one-loop effective 
action. I. e., it will be obtained as a self-consistent solution of the one-loop held equation: 

The problem is, that no analytic representation of the effective action in the one-loop approximation is known. In 
pT| a numerical method to calculate the functional derivative of an effective action has been presented. It is based 
upon a computation scheme to determine full one-loop Green’s functions 0-Q. In this paper this technique is 
applied to determine a bounce-solution from a held equation where the full space-time-dependence of the non-local 
terms of the quantum corrections to one-loop order is taken into account. As the numerical method to compute the 
functional derivative needs the background held as a numerically known input quantity the procedure we use to obtain 
the background field is iterative. It finally results in a self-consistent fixed point of the iteration condition, i. e. a 
solution of the one-loop field equation (||). 

A one-loop Green’s function is a divergent quantity that has to be renormalized. An important issue of the presented 
method is the strict separation of the divergent parts, which have to be handled analytically, and the convergent part, 
which may be computed numerically. As no expansion in any quantity is needed a conceptual advantage of this 
method in contrast to approximation schemes like the derivative or the heat kernel expansion is that the unphysical 
infrared divergencies these approximations suffer from do not appear. 

Among the various models for which self-consistent solutions are interesting this paper deals with the electroweak 
phase transition. Recent lattice calculations show that the perturbation theory is not reliable for Higgs boson masses 
above 60GeV. The phase transition seems to disappear at a critical Higgs boson mass which is in the range between 
66GeV and 80GeV and the experimental lower limit of the Higgs boson mass already is beyond 66GeV. Hence, 
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the generation of the baryon asymmetry at the phase transition of the standard model is ruled out already, and 
the interest has shifted to extensions of the standard model. Nevertheless, the investigation of the first order phase 
transition at small Higgs boson masses is still of interest and is at present actively studied on the lattice |2^]. Some 
aspects of the transition, like the determination of the surface tension and the latent heat, are still in a state of 
development. Hence, we have the possibility to develop new calculational methods in a well known model that are 
straightforward to extend to more complicated models. 

The plan of this paper is as follows: In the next section the model will be specified and the relevant fluctuation 
operators are given. The next section deals with the renormalization of the effective action that is improved by 
resummation. In section IV the formulas for the numerical procedure how to calculate the functional derivative of an 
effective action and the effective action itself are derived. Finally the results are presented and discussed. 


II. BASIC EQUATIONS 

The Lagrangian of the SU(2)-Higgs model with gauge fields and fermions in the limit of vanishing electroweak 
mixing angle is given by 


^ — -^gauge -^Higgs .^ferin j 


where 

/:gauge = , (7) 

/:Higgs = ^(D^<I>)t(D'^<I>)-t/(<i>^$), (8) 

/:ferm = ^ ^ , (9) 

/ //' 

it;, = d^Al - , (10) 

, ( 11 ) 

C7($t$) = ^ ^ (12) 


Here W;, denotes the field strength tensor of the gauge fields, $ the Higgs-doublet, t“ the Pauli matrices, the 

f f' 

fermion fields of the quarks, uq the vacuum expectation value of the Higgs field at temperature T = 0 and gi/ the 
Yukawa-couplings of the fermions to the Higgs field. The sum / runs over the different flavors and colors of the 
quarks. 

The critical electroweak bubble is a pure real Higgs field configuration. On the classical level it is 


$(x)ci = 



= 0 , 5 '(a:)ci = 0 . 


(13) 


The configuration is taken to be stationary, i. e. time independent. As explained in the introduction the bubble profile 
calculated from the functional derivative of the classical action (^) does not exist because the tree-level potential U{4>) 
only has a single minimum. Hence, instead of the classical action one has to use the effective action. For the model 
under consideration this effective action is given by: 


T[<P, T] = + S'eff.gauge[<^, T] + S-eff.ferm [<(>, T] + Higgs T] . (14) 

This functional cannot be computed exactly. Usually it is calculated in the semiclassical approximation to one-loop 
order. In order to calculate the one-loop approximation one needs to know the fluctuation operators. They are found 
by expanding the Lagrangian (|^ around the classical field configuration ( p^ ) to second order in the small fluctuations: 

$(x) = [(/)(a;) -f h{x) + *r“(^“(a:)] , 

=<(a;) , (15) 

4'(a:) = ipix) . 
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Here h{x), a“(x) and i^ix) denote the fluctuations of the isoscalar- and isovector-part of the Higgs held, the 

gauge field and the fermions respectively. The second order contribution to the Lagrangian is given by: 

f If 

+^d^hd^h + (16) 

-i (3A(()2 - Av^) ^ £“(/?“ , 

One has to add a gauge-fixing term and the contributions of the Faddeev-Popov-ghosts to the Lagrangian. Here the 
‘t Hooft-Feynman-background gauge is used: 

with jra ^ ^ (^ 7 ) 

/Iff = ^7^“ »?“ • (18) 

From = (£ + £gf + £fp)^^^ the fluctuation operators are read of. They are given in Euclidean metric as the 
bounce solution we want to obtain is a solution of an Euclidean field equation. There are three contributions of the 
bosonic sector to the effective action. The first is due to the isoscalar fluctuations of the Higgs field, the second is 
due to the coupled channel of the gauge-fields and the isovector components of the Higgs field (would-be-Goldstone 
bosons) and the last is due to the time component of the gauge fields and the Faddeev-Popov-ghosts. Each of them 
can be written as: 




Ci In Ji = a In det 


( —92 -|- \ 

\-d^ + u,{o)) 


In order to shorten the notation the quantities 
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/ 

0 
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0 
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0 

0 


0 

VO 
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0 

m2(())) 

(C^ 

<Y. 




+ - Up) 

\ 

TKW , 


/ 0 0 0 1 \ 
0 0 0 0 I 
0 0 0 0 
V 1 0 0 0 / 


/O 0 0 0 \ 

0 0 0 1 3 
0 0 0 0 ’8 
Vo 1 0 0/ 


/O 0 0 0 \ 
0 0 0 0 I 
0 0 0 1 
Vo 0 1 0/ 


(19) 


( 20 ) 

( 21 ) 


( 22 ) 

(23) 


are introduced and one uses that the background field is time independent. 

The contribution of the top-quark is the dominant one in the fermionic sector. Therefore, nothing but its contri¬ 
bution is taken into account. It is given by 


= -3 In det , (24) 

where the factor three is due to the quark color and gt is the strength of the Yukawa-coupling of the top-quark. As 
explained in the introduction only the leading temperature dependence of the top-quark contribution, which is part 
of the finite temperature effective potential, will be taken into account. Hence, this calculation is not self-consistent 
with respect to the fermionic sector. 
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The first step is to extract the contributions of the different fluctuating fields to the high temperature approximation 
of the one-loop effective potential [| Based upon the fluctuation operators given above the sum of the tree 


level potential C/(0) and these one-loop corrections can be written as: 


with 

D 

B 

rp‘2 

-^0 

E 

Xt 


- ET4>^ + 

2D{T^ - T^) 

(3m^ -I- + 2m() 

(2m^ + [-f- + ml + ^-m^ - 

ml — 8 vqB 
4D 


327r 

A- 






2^4 In^ 


-h 


m 


H 


I'H 


4 , 




m 


w 


In 


TO- 


W 


asT^ 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 


where Inos = 21 n 47 r — 27 ~ 3.91, InoF = 21 n 7 r — 27 ~ 1.14, vq = 246GeV, mw = ^gvQ and ©w = 0. The form of 
the constants B, D, E and differs from those in [y, because of the contributions of the Higgs boson and would- 
be-Goldstone boson fluctuation taken into account here. The renormalization has been performed in the minimum of 
the classical potential at T = 0. 

The given representations of the effective action are quite formal. They are divergent and have to be renormalized. 
Furthermore, they have to be evaluated for a system at finite temperature. This is done in the next section 


III. RENORMALIZATION OF THE RESUMMED EFFECTIVE ACTION 

The vacuum structure is changed due to the contributions of the different fluctuating helds to the one-loop effective 
potential. The tree-level mass and the tree-level self-coupling of the Higgs held are modihed due to effects of hnite 
temperature. They will be modihed again by effects of higher order. Dependent on the temperature these corrections 
are not negligible and a pure one-loop calculation is not reliable . 

In order to take into account at least the dominant contributions of the higher order effects one has to consider 
the sum of one-loop-Daisy-graphs, replacing the expressions due to the inserted tadpole-graphs by their leading 
temperature dependence which is proportional to . Also the vertices are replaced by those which take into account 
the temperature dependence of the coupling. 

In practice the summation of Daisy-graphs is performed by replacin g th e tree-level-potential U{(jj) in the Lagrangian 
by the high temperature potential I4t given in the preceding section ||3l|-M■ Then the fluctuation operators and the 
contributions to the effective action are determined. The effective action based upon these fluctuation operators is 
improved by resummation. It can be written as 




n— — oo 


where lJi{4>,T) is given by (^ff. with TO^((/)) and m^{(j)) replaced by 

o2 

m\{(j), T) = — &ET(j) + ml , m^^{(j), T) = -^4’’^ + — 8ET(j) + ml . 

Expanding the effective action in a series of one-loop Feynman-graphs it can be rewritten as 

E Tr{[-A-t-:.2^Z^,(0,r)]“'[W,(<^,r)-W,(0,T)]} 

k—1 n— — oo 


(32) 


(33) 


(34) 
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(35) 


+ 00 

=E-5.‘5^ 

k=l 

Here [—A + -\- Ui{Q,T)]~^ is a formal representation of the Green’s function associated to that operator. The 

superscript (fc) denotes the order of the corresponding graph (see fig. 1). 

In analogy also the improved one-loop effective potential T) exists. It has the following relation to the usual 

one-loop effective action I4((/', T): 


r d^D 

T) = c,T tr ^ ^ In [p^ + ^1+ Ui{(t>, T)] 

-|-CXD p ^2 

= c,T tr ^ / -7^ In [p^ + vl+ 

T) 

+ vl+Ui{4>) _ 

= I4((^,r) + H,,ring(<(>,T) . 


-l-CiT tr ^ 


■ d^p 


In 


(36) 

(37) 


(38) 


The divergencies in Vi are temperature independent. They are cancelled by the familiar T = 0-counter terms. 
However, also V^^ring contains divergencies which furthermore are temperature dependent. Hence, the substitution of 
U{(j)) by Vht(<)>) accounts for finite and infinite contributions of higher orders. The latter are removed by the use of 
thermal counter terms They are not fixed by renormalization conditions but are chosen in such a way that 

they cancel the two divergent orders of Hi_ring(<)>, T) exactly. These two orders (fc = 1 and k = 2) can be read of from: 


+ 00 


+cx> 


V^z,ring(</>, T) = ^ H^tg = C.T tr ^ 


d^p ^ (-l)fe+i 

Jc 




(27r) 




p^ + vl+U^{Q,T)\ 


p^ + vl+U^{S),T)\ 


(39) 


f 2) 

However, the high temperature contribution of Voting; i- contribution of the static Matsubara mode, is finite 

and is not subtracted. 

For the numerical application the resummed one-loop effective action is separated into a finite part that is evaluated 
numerically and an infinite contribution that has to be renormalized. The latter is given by the first two orders of the 
series (^). In performing the renormalization by adding the counter terms one finds contributions already taken into 
account in the high temperature potential t4t(</', T); we denote these terms as I^,ht- These terms must be subtracted 
in order to avoid double counting; we call the resulting difference 


A c(H-2) _ ay 


:(i) 


s?l-s, 


( 2 ) 

i,R,ht 


-P 


I d"* {Vl 


(1) 

ring 




( 2 ) 


-K 


( 2 ) 


i,ring,ht 


K 


ht , 


s^, 


(40) 


( 2 ) 

Here we make use of that the high temperature approximation of the second order contribution, is finite. 

Therefore, it is treated together with the finite part, and is evaluated numerically: 


C _ 0(2) I 0(3) . rr(3) 

‘i’r.R.num — Cij -|- Ci- j, 


= 5, 


( 2 ) 

i,R,ht 


s} 


(3) 


z,R,ns 


(41) 


Here a superscript (k) denotes the sum of contributions beginning with the fc-th order. The index ‘ns’ denotes the 
contributions of the non-static Matsubara modes. It is known to be approximated sufficiently by a gradient expansion 
It is found to be negligible small compared to the contribution due to the static mode. Therefore it is neglected 
( 2 ) 

and only is evaluated numerically. This is postponed to the next section. 

Using 


P ^ f Zz _ 1 ^ [ d'^p 1 

Z' J (27r)3 p'^ + Vn J (27r)^ p^ -|- rrP 


(42) 
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12 47r 167r2 asT"^ 

d'^p 1 

(27r)^ (p2 -(- rri^Y 


+ 5i{T,m^) , 


( 43 ) 


In- 


87rm IGtt^ obT^ 


+ (52 (T, m^) 


and neglecting the terms with (5i and (52, which are small at the temperature we are working at, the analytic part 
ren found to be: 


A4';e?=/3 


+ 




327r^ 

g^T 
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327r2 

5^ 
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^(ln-^ + 2) - m^l - 


2 l^''aBT2 
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1 ?a-}' 


(44) 


IV. DESCRIPTION OF THE NUMERICAL METHOD 
A. The functional derivative of the effective action 

In the previous section the calculation of the finite contribution to the resummed effective action due to the different 

fluctuating fields was postponed. Combining them to its functional derivative with respect to the background 

field (j) can be represented as 


(5.S*R,ht 

6(j}{x) 


= tr 


SU{cl>,T)^,^ ; 

C —Xj 


= tr 


5(j){x 
(57W^((/>,r) 
5<j){x) 


G{x, x) 


— giv[c CVG(a;, x) 


(45) 


where W, T) and C are the (6 x 6) generalisations of the matrices defined in ([2T|)ff. at finite temperature. They 

are given explicitly in appendix |^. c contains the factors of degeneracy for the different fluctuations. Its non-zero 
components are given by 


cii = C 22 = C33 = C44 = —C66 = 2 > and 


C 55 = 2 . 


The Green’s function G{x,x') is the solution of 


-A+G((/),T) 


G(a;, x') = 1 • (5^^^ {x — x') . 


(46) 


(47) 


A method for the numerical computation of such a one-loop Green’s function was derived in p2|-p4|. In order to 
determine G{x, x) it is decomposed into its partial wave contributions. The details of this calculation are given in 
appendix It turns out that by the decomposition one channel of the coupled system can be combined with the 
Faddeev-Popov channel so that the system reduces to a (5 x 5)-system. Extracting the order (2) from ( |4^ ) it takes 
the form: 


i;c(2) - 1-00 

S^{x) - ^ 


21 + 1 \ 3g^ 


47r 


(^iil + 522 !;) + 2 ( 


+3{XTcl)-ET)gi\]]-^ct>gg], 


2 E 0 ^ i ^ 13 .) + ^9i3,i j + <^0 ( 523 .) + ^9^,1 


21 + 1 


'( 1 ) 


( 48 ) 


hi) 


'( 1 ) 


( 1 ) 
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Here the radial Green’s functions gkm,i{'r,r) are defined by: 


Mnkgkm,iir,r') = -^5{r - r') 

with M := M° + 1/ , 

^2 


M" = + ml 


m„ = (0 , 0 , friH , fhH , 0)„ , 

= (^ + 1 , ^ ~ 1 , I , I , l)n ■ 


The non-zero components of the potential V (</>) are 


Vn = = V"55 = , V33 = ( ^ + At ) </>'^(r) - 3 ET<f>(r) 


V44 = 3AT</>^(r)-6ET</^(r) , V43 = V31 =- 

and F 23 = VJi 2 = (r) ■ 


1 + 1 
21 + 1 


g(t)'{r) 


(49) 


(50) 


This coincides with ||^ where it was derived by reduction of the sphaleron system. The radial Green’s function is 
given by [p§ 


gkm,i{r,r') = K /fet(?’)/“^;(r')0(r-r') + {r)f^'^i{r')Q{r'-r) 

where f^f and are the linear independent solutions of the homogeneous equation 


2 d Uln + l) 2 


r dr 


- Snk - Vnk{4>,T) 


4T W = 0- 


(51) 


(52) 


The mode functions {fki ) regular (singular) at infinity. In order to calculate the mode functions numerically 
their asymptotic behaviour at infinity is separated from the rest by the following ansatz 


The functions bf^ (z) satisfy the equation 


with Z = KT . 


/ 92 2 d Uln + 1 ) 

\dz^^ zdz z 2 




One obtains 


b+{z) = z“('+2 )(2Z -h 1)!! 

{2l-l){2l + l) 

K(^) = - 

15 2 




(2;-h3)!! ’ 
z^-\2l + l) 
{21 + 1)11 ’ 
KM) = M) , 


KK) = 1 • 


(53) 


(54) 


(55) 


In the massive channels we obtain the modified Bessel functions, while in the massless ones the solutions reduce to 
powers of z. The fifth channel only contributes in the s-wave. Therefore, its Bessel function is given only for I = 0. 
The normalization constants and k = ttih are chosen in such a way that the Wronskians of the different channels 
become equal: 


W 


KM’K(z) 


(56) 


















The differential equation of the rest mode functions reduces to: 


r ^ 

^ 1 

5 ] 

dr2 


dr 


Kiir)=Vnk + 




(57) 


This equation can be used for generating the functions order by order in the potential V. Introducing the 


contribution of order k in the potential as and defining 


_ +00 

_ 


-EC 

j=k 


±U) 


as in 0 , the order (1) of the Green’s function is given by: 


(58) 


B. Evaluation of the effective action 

In order to calculate the nucleation rate 7 the effective action in the one-loop approximation is needed. A method 
to calculate its functional derivative with respect to the background field has been given in the previous subsection. 
For the purpose of renormalization the divergent parts have been separated from the rest and have been handled 
analytically. The same separation can be used here in order to evaluate the effective action itself. The analytic part 

is given by (0)ff. The numerical contributions are determined in the following. 

For the numerical computation of the effective action the fluctuation determinant J is decomposed into partial 
waves 


lnJ = ^(2Z-bl)lnJ, , 

i 


with 


det M 
det MO ■ 


(59) 


Here M and M*^ are the fluctuation operators defined in (p^. We need to know two (n x n) matrices fi(!z, r) and 
f°(z/, r) which contain the n linear independent solutions and of the system of coupled differential 

equations 


+ and + = (60) 

The boundary conditions of these functions are chosen to be regular at r = 0. The lower index denotes the n 
components while the different solutions are labelled by the Greek upper index. These solutions are normalized such 
that 


hhn/i(i/, r) (/;°(i/, r)) 

On these conditions the statement of the theorem is 

, ^ det (M -t- I/O 

det (MO -b 1 / 2 ) ^ ™ 


= 1 . 


det f{i 


det /°(j/, r) 


(61) 


(62) 


The theorem has been applied for computing the one-loop effective action of a single scalar field on a bubble background 
[p0| , and of a fermion system on a similar background 191 previously. Furthermore it has been used to calculate the 
bosonic fluctuation determinant of the critical bubble |17|. It was found to yield very precise results, in addition to 
providing a very fast computational method. 


9 















If the theorem is applied at i/ = 0 the fluctuation determinant J = J{Q) is obtained. The consideration of finite 
values of u is necessary in the discussion of zero modes. For simplicity the following formulas are given only for v = Q. 
In analogy to the mode functions f^{r) also the functions fi{r) can be decomposed as: 

fkA^) = , (63) 


with 6;^ {kt) given by (|^). The consequence is that the separated functions satisfy the same type of differential 

equation as the rest mode functions /i^y(r). Only the boundary conditions differ, namely h^;(0) = 0 and h'^ i{oo) = 
const. The two functions are related by: 


hi {r) 


1 + hi{r) 


1 + hi{oo) 


- 1 . 


(64) 


After determining fi(r) one obtains the desired order of the fluctuation determinant. 

InjT”/^^ = lim |ln det[l +/i;(r)] — tr| . (65) 

As the determinant J is a product of three individual determinants due to the coupled and uncoupled channels one 
has to take into account the different factors of degeneracy for each of them (see (|4^)). 


C. Self-consistent determination of the bubble profile 

1. Computing the start profile for the iteration 


In order to solve the one-loop held equation iteratively we have to choose a suitable separation of the equation into 
a left-hand and right-hand side. Then setting the right-hand side to zero results in a first approximate solution of the 
equation. This is taken as the start profile for the iteration. 

In particular we choose the following form of the held equation 
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3 Xi71hT 3A 
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Stt 

3mHT 
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_ (ffA 
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Utt 

2 
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SitihT 3{jn\j - rri^) 


asT^ 
2 


Stt 


327r2 


3mu , m 


-k At U - 3ET(I) 


w 


647r2 asT^ 


+ 2 - 




( 66 ) 


(67) 


I. e., in addition to I4t the contributions of the isoscalar and isovector part of the Higgs field to the resummed effective 
potential are shifted to the left-hand side. For a suitable chosen temperature the potential Vfi has a secondary minimum 
at a nonzero value of We denote this minimum as v(T). It is used to scale the dimensional quantities as follows: 


(j) = v(T)(j) , Xfj, = J, , and T = gv{T)T . 
gv{T) 

All masses are scaled with gv{T). The inhomogeneity T((/)) consists of (Q) and (|4^ ) 
shifted to the left-hand side. 


( 68 ) 


reduced by those terms that were 
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(69) 


n-f ^ 1 / (1) , W \ , f 9^ , n\ M o J^rr\ ^ 2/ + 1 ( 1 ) 

A^) = [9n,i + 9^1 j + ( + 2 At<(> -3ET n -^933 

I 7—n V / 7—n 


Z=0 


S^r ^27 + 1 

^ 27r 

1=0 


'( 1 ) 


,(1) 


I 9i3,1 + 9i3,1 ) ■*" '^0 ( 923,1 + „923,1 


'(1) 


,(1) 


+ 00 


-3T (At(I)-ET)J2 

3\TT(j) 


27 + 1 ^ , 3g^T^ —) 

- 9h,i + 


47r 


47r 

3T ^ 

Stt 


1=0 

(3XtA + ^h) ^ 

,2 


167r 

39^T^^ 3g^ 


IGtt 


327r2 


A())ln 


m 


w 


asT^ 


+ At U - 3£;T 


ff 


+ Xt ] <p — 3ET(j) + fnjj I — niH 


For the purpose of comparison the solution (j)c\ of the equation 

— A(/)ci + A 


= 0 


0 = 0cl 


(70) 


is also calculated. It is denoted as the classical profile. 

In order to avoid large changes between two subsequent profiles especially at the beginning of the iteration we use 
a relaxation, i. e. the field equation (p6|) is written as 




= (l-e)Z((/)(”-i))+£!((/)(”)) . 


(71) 


By a suitable chosen e the convergence of the iterations of the profile can be improved. 


2. Treatment of zero modes and unstable mode 


For the isoscalar fluctuations there are two contributions to the effective action and its functional derivative that 
have to be treated individually: namely the instable mode in the partial wave 7 = 0 and the translation zero mode in 
the partial wave 7 = 1. 

A bound state occurs in the s-wave, i. e. a state with a negative eigenvalue. As this is the only negative eigenvalue the 
total determinant is also negative and the one-loop contribution to the effective action is complex. This corresponds 
to the fact that the state under consideration, the critical bubble, is unstable. 

In the p-wave there exists a mode with the eigenvalue zero. It corresponds to the fact that the critical bubble 
breaks translational invariance. Because of this mode the determinant itself is zero and the effective action diverges. 

For the calculation of the nucleation rate instead of the fluctuation determinant with these two properties we need 
one where the negative eigenvalue has been replaced by its modulus and where the three eigenvalues zero have been 
removed. These changes to the determinant are indicated by two primes: 


+ 00 

A= n 


n= —00 


det" [-A + I/2 +m^((/),r)] 
det [-A + jz^ +m2(0,r)] 


(72) 


As three eigenvalues zero have been removed from the determinant of the isoscalar fluctuations the rate gets its 
dimension (energy)^, i. e. (volume)”^. The further dimension (time)“^ is due to the eigenvalue co- of the unstable 
mode. 

The numerical method presented in the previous sections does not explicitly make use of the eigenvalues. Therefore, 
we have to show how these two contributions are removed from the effective action ~ this has also been discussed in 
I p^ - and its functional derivative. 

In case of the unstable mode this is performed easily. One computes the particular partial wave contributions 
according to (^). In the partial wave 7 = 0 only the sign of the determinant has to be inverted before taking the 
logarithm. For the functional derivative the sign is irrelevant because due to 


S 

6(j){x) 


InJ 


S 

6(j){x) 


InJI^^ 


E 


uXi 6(l){x) 


(73) 
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the sign of each eigenvalue enters the equation twice. 

The removal of the translational zero mode is more complicated. An alternative representation for the Green’s 
function 544 _i(r, r') compared to (^l|) can be deduced from the eigenvalues and eigenfunctions of 




r dr 


i^a{r) = ujl'ipc.ir) 


(74) 


It is given by: 


ff44.i(r, r') 


V’a(OC(^') 

rv ^ 


(75) 


Obviously this Green’s function is divergent when the spectrum of the operator contains the eigenvalue zero. Defining 
a new Green’s function as 


544,i(r,r',i/) = ^ 

Ct 


^a(r)C(r') 
+ UJa 


(76) 


the divergence manifests itself in the limit z/ —> 0 in the form of a 1/z^^-pole. As the zero mode is the only mode 
which contributes to this pole its contribution can be removed by defining the Green’s function at = 0 without the 
zero mode in the following way: 


gi 4 :,i{r,r', 0 ) := lim 

i/^O 


544,z(r, r',z/) 


7/2 


(77) 


In the special case under consideration, where the Green’s function only has to be calculated at r = r', the practical 
way to subtract the zero mode in the p-wave is to compute the function 544 p(r, r, z^) for sufficient many different 
values of ly. Then in the function 


544.1 (r, r, v) 


^ + B(r) + C(r)z.2 


(78) 


the parameters A, B and C are adjusted by a fit with respect to v for each value of r. The subtraction of the zero 
mode is then done by setting 544,1 (r, r, 0) := B{r). 

The computation of the function 544^1 (r, r, v) is done in the same way as described in section tV. The only difference 
is, in the definition of the radial part of the Green’s function ( ^ ) and in all consecutive equations a has to be 
added to the partial wave operator M°. Hence, the ‘mass’ k of the isoscalar channel changes to zt = yuh^G-zA. 

The one-loop contribution of the isoscalar fluctuation to the fluctuation determinant is determined according to 
(p^ in each partial wave. In the p-wave the following expression has to be evaluated: 


lim In 

r—*oo 


l + hli{r) 


(79) 


The numerical investigation shows that hf 4 ( 00 ) = — 1 and the logarithm diverges due to the zero mode in this partial 
wave. In order to subtract the mode the computational method to determine the fluctuation determinant is extended 
in the same way as before in case of the Green’s function. I. e., the theorem mentioned above is applied at finite v. 
In order to subtract the zero mode the determinant Jh,i{v) has to be replaced by |p7|,p0| 


lim lim 

V —>-0 r —>^00 



(80) 


In practice h\ ^(z/, 00 ) is computed for some sufficient small values of v and is used to fit the constants A, B and C in 


1 -I- h\ i{v, 00 ) = Av'^ + Bv‘^ -\- Cv^ . (81) 

With these techniques the unstable mode and the translation mode are removed from the effective action as well as 
from its functional derivative. Concerning the translation mode a further remark is necessary. An exact zero mode 
only exists when the profile (j) is determined from the classical field equation ( [^ . The one-loop equation differs from 
the classical one due to the non-vanishing right-hand side of the equation. As a consequence one finds a bound state 
in the spectrum in the p-wave below the continuum threshold but its eigenvalue only is near to zero but it is not 
exact zero. Nevertheless this mode is subtracted. 

In the coupled system of gauge fields and would-be-Goldstone bosons a zero mode exists in the s-wave. For the 
treatment of this mode we adpopt the prescription given in |l7|| . 
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3. Some details of the numerical evaluation 


In order to determine the profile for one step of the iteration a boundary problem has to be solved. At r = 0 the 
slope of the profile has to vanish and at r = oo the function has to tend to zero exponentially. The profile is obtained 
using a shooting method: The value of the profile at r = 0 is varied until the profile gets the correct behavior at 
infinity. We typically use 1000 steps on a scale from r = 0 to r^ax ~ l/w// where the latter typically is below 200 in 
the units gv we are working in. The numerical technique used to integrate the differential equation is the Nystrom 
method. 

For small values of r the profile can be represented as (j){r) ss oq + a 2 r^ + a^r'^. Hence, for a given profile the 


parameters can be obtained using a simple fit. Then for the rest mode functions which are calculated from the 
differential equations (57) the behavior at small r can be evaluated analytically using a power series expansion. Using 
this expansion the functions have the correct starting properties for the following numerical integration. Hence, we 
can avoid numerical problems due to the divergent behavior of the Bessel functions in this region. 

The radial Green’s functions are composed from the rest mode functions according to (^) and then the partial wave 
summation (M) is performed. This sum is executed explicitly up to /max = 15- At large / the product (2/ + r) 

behaves as + b{r)/P + c{r)/P, no matter which channel i one considers. Using a fit the coefficients a, b and 

c can be obtained for each value of r and then the value of the sum from /max to infinity can be estimated. As a 
cross-check we start this estimate at / = 10 and it is tested whether it approaches a constant when / reaches /max- 

In order to check the precision various cross-checks have been implemented. The equality 


:= rfiHr'^W 




n^l 


= 5^ 

ru 


(82) 


is satisfied numerically with a relative deviation of the order 0(10 ®). Especially the Wronskian between f~^{r) and 
/(r) satisfies: 


KT^W 


r+(r),/^(r)] =5i + hii^) 

= S^ + /i“+(0) . 


(83) 

(84) 


The relative deviations from these identities are of the order 0(10 for (^) and 0(10 ®) for (| 8 ^). The former being 
not as precise as the latter is due to the bad convergence of the functions at large r. In the massless channels they 
only achieve their asymptotic values as l/r7 Hence, their value at r = oo only can be estimated from a fit. On the 
other hand the constancy of the product of with the Wronskian is satisfied excellently. Therefore we do not make 
use of this fitted quantity in the further process. 


V. RESULTS AND DISCUSSION 

The numerical results are given in table 1 and in the figures 2 to 7. The Higgs boson masses used here are 30GeV 
and 40GeV. The mass of the top-quark is taken to be to* = 170GeV and for the gauge boson mass we use a value 
of mw = 80.2GeV. The gauge coupling g is taken as g = .651616. The vacuum expectation value of the Higgs field 
is then fixed to be vq = 246GeV. The relevant temperature interval is determined from the approximation of the 
resummed effective potential Vr. It is found to be smaller than the interval determined from the high temperature 
approximation Uht of the effective potential. The values are given in table 1. Within this interval we choose equidistant 
steps for the temperature. We leave out the range of temperatures close to the roll-over temperature Tq, because 
these correspond to small thick wall bubbles which are of minor cosmological interest. 

We give the results for the logarithm of the nucleation rate 7 , i. e. the rate including the one-loop corrections and 
prefactors as defined in (|^). It is evaluated using two different bubble profiles. First we use the self-consistent profile, 
i. e. the solution of (^). The results have to be compared to that ones we obtain from using the profile which is 
the solution of The results are displayed in the figures 2 and 3 where we give the logarithm of the nucleation 
rate versus the temperature. For the results based on the self-consistent profile we also present a function that fits 
the explicitly calculated points very well. The function has the form 

'’(^1 = -»( 7 #^)* ■ 

The results for a and b depend on the specific parameter set. We find a = 3.48 • 10“^ and b = 1.350 for rriH = 30GeV 
and a = 2.38 • 10“^ and b = 1.362 for ttih = 40GeV. 
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At ttih = 30GeV the self-consistently calculated transition rate is smaller than the rate obtained from the classical 
profile, while at mn = 40GeV it is just the other way around. This is due to the relation between the temperature 
intervals. While in both cases the interval suitable for finding the self-consistent solution is located within the interval 
derived from the high temperature potential it is shifted towards higher temperatures with increasing Higgs boson 
mass (see table 1). For comparison we map the temperature intervals [Tj^, T^] and onto the interval [0,1] 

by introducing the dimensionless quantity ( = {T — To)/{Tc — Tq). The results of both parameter sets are given in 
figure 4. For both parameter sets the self-consistently obtained rate is increased compared to the ordinary one. Hence, 
the effect of using the self-consistent profile is an enhancement of bubble nucleation. 

Based on the results for the nucleation rate we can investigate some aspects of the cosmological electroweak phase 
transition. As the universe expands the temperature T decreases. Once the critical temperature is reached critical 
bubbles can start to nucleate and grow. The corresponding time is obtained using the relation between time and 
temperature 

t ~ m-pi, = 1.2 • 10^®GeV . ( 86 ) 

Hence, the nucleation rate 7 can be considered as a function of time as well. For a given velocity v of the bubble wall 
the number density of bubbles of size r at the time t is given by 




Then the fraction of space that is still in the symmetric phase is obtained as 


Ps{t) = exp 


I 

-y J dti 


(87) 


( 88 ) 


This quantity is used to define implicitly the time te when the nucleation is completed by Ps(te) = 1/e. As it marks 
the end of the phase transition it is a quantity that is easier accessible in lattice simulations than other characteristics 
of the transition. The fraction of space that is already in the asymmetric phase is simply Pa{t) = ^ — Ps{t)- The total 
number of droplets per unit volume, N(t), also depends on the fraction Ps(t) because new bubbles can only nucleate 
in a region that is still in the symmetric phase: 


t 

N{t) = Jdtij{ti)ps{ti) . 

tc 


(89) 


This quantity is important for the discussion of structure formation. The average bubble radius, R{t)^ is obtained as: 


t 

J dti v{t - ti)-/{ti)pa{ti) . 


(90) 


The quantities Pa{t), N{t) and R{t) are evaluated for the parameter set uih = 30GeV, TOj = 170GeV for different 
values of the bubble wall velocity v. The results are displayed in figure 5 to 7. The completion time depends only 
weakly on the bubble wall velocity and is roughly te — tc ~ 0.15 • 10“^^sec (see figure 5). This can be compared 
to a previous calculation of this quantity [Q, where a value of roughly 3 ■ 10“^^sec is given. One sees the effect of 
enhancing the rate. The time scale on which the nucleation proceeds is decreased by an order of magnitude. Also 
for the number of bubbles per volume after the completion time the results can be compared. For three values of 
the bubble wall velocity, namely c, c/10 and c/1000 the corresponding values for ln(A^/cm^) can be read of from 
figure 6: they are 28.1, 34.7 and 47.8. For the same velocities in [Q the values 9.3, 12.3 and 18.0 are given; hence, 
the number of bubbles per volume is increased substantially. The average bubble radius R{t) for the three velocities 
under consideration comes out to be the same as in |38| and is plotted in figure 7. As the initial radius of the critical 
bubble is of the order of Ifm it is neglected throughout this calculation. 

The presented method to calculate self-consistently the solution of a one-loop field equation has the following 
essential features: 


• The complete one-loop determinant and its functional derivative is evaluated. No expansion in any quantity 
is needed. In view of the fact that the derivative and heat kernel expansion generate unphysical infrared 
divergencies order by order this is an important conceptual advantage of the method presented here. 
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• Coupled with an iterative procedure the numerical computation of the functional derivative of an effective action 
allows the determination of a self-consistent solution of a full one-loop field equation. 

• The divergent and convergent contributions to the effective action can be treated separately. While the conver¬ 
gent part is calculated numerically the divergent part is evaluated analytically. So the choice of the regularization 
and renormalization prescription is free. 

• The numerical computation of the functional derivative of an effective action only requires numerical integration 
of a few differential equations per partial wave. This can be done with high precision and is not computer time 
consuming. 

In conclusion we find that taking into account the radiative corrections already in the determination of the bounce 
solution has the following effects: 

• The width of the temperature interval within that a first order phase transition is possible is decreased and with 
increasing the Higgs boson mass it is shifted towards higher temperatures. 

• The transition rate is increased substantially relative to the rate obtained from the standard calculation. 

• The time scale within that the transition is completed is reduced. The number of bubbles per volume is increased. 
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APPENDIX A: PARTIAL WAVE DECOMPOSITION OF THE GREEN’S FUNCTION 


In the partial wave decomposition of the Green’s function was done by reduction of the sphaleron system. 
An alternative derivation will be given here. The Green’s function G{x, x') introduced in (^) satisfies the matrix 
equation 


X>) = 1- d^^\x - x') . 

The non-zero components of T) and C are given by 

Ml^ = MI 2 = Mh = Mh = mli^) , Ml = m2 (</>, T) , = mU<p, T) , 

Cl4 = Cll = 1 5 C24 = C42 = 1 ) C34 = C43 = 1 ■ 


(Al) 


(A2) 


For the partial wave expansion of the solution 'I'(a;) of the corresponding homogenous problem we use the following 
ansatz: 


+ 00 

1=0 


ll)l{x) = 


+l 


m— — l 


(ip\{x) iflix) ip\{x) 0 0 0 \ 


( hir) \ 

lpI{x) ipl{x) ipl{x) 0 0 0 


f2ir) 

lpIIx) (filix) 0 0 0 


foir) 

0 0 0 ip\{x) 0 0 


fir) 

0 0 0 0 ^l{x) 0 


hir) 

\ 0 0 0 0 0 ifilix) ) 

\ U{r) / 


ip]{x) := [Yi 1+1 m{x)]^ , ^]{x) := [Y 1 ;_1 m{x)]^ 

ipi{x) ■= ipl{x) = ipl{x) = Yi^{x). 


iplx) :=[Yii m{x)]^ 


(A3) 

(A4) 

(AS) 

(A6) 


Here x := x/r means the unit vector in x-direction. To shorten the notation the I- and m-dependence of the 
functions tpl and the ^-dependence of the functions fi is not shown explicitly in the following. The vector spherical 
harmonics Y (see e. g. [p^) are related to the common surface spherical harmonics: 
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Yi 


/+1 m — 


^ I 1 — 1 m — 


1 1 

Cl 2 / “h 1 

1 1 


Cq 2^ “h 1 

YI I ra = — 7 ===ILyi 


[ixYiTn + rVYim^ , 


\/l{l + 1 ) 


I 


(A7) 
(A 8 ) 
(A9) 

(AlO) 
(All) 

For the partial wave I — 0 one has to consider the following peculiarity: neither the function IFo -io(®) nor IFo o o(®) 
exists, because the normalizing prefactor diverges. Scaling / 2 (r) with cq and fsir) with ^l{l + 1), the products 
Yi i-i mf 2 {r) and Yi i mfsir) are well defined also at Z = 0; they yield zero. I. e., at I = 0 the components f 2 {r) and 
fsir) vanish identically. Hence, the coupled (4 x 4)-system reduces to a coupled (2 x 2)-system. 

In the following the index I of the functions ip is suppressed. Using the index notation the ansatz ( |A4D can be 
formulated as 


Co := 


21 + 1 


Cl := 


1 + 1 
21 + 1 


and with 


T^nk ■ = 


i’n = +nfi , 

2 d 1 


^ dr"^ r dr r 

the homogeneous equation for dt takes the form: 


- -TL^ ]Sr.k-M 


2 

nk 


Using 

we get 


'Dnk'fpk = g(p'Cnk—^k ■ 


TL‘^(p’^ = li{li + l)(p'!^ with k = {I + 1,1 - l,l,l,l,l)^ 


TdnkP^k — +k 


2_5 kik + 1 ) 
dr^ r dr 


^nk AI 


nk 


h = g^'Cnk-^A • 


(A12) 
(A13) 

(A14) 
(A15) 

(A16) 

(A17) 

(A18) 

(A19) 
(A20) 

Therefore, as the vector spherical harmonics are linear independent, ( A16|) only can be satisfied for fk being a solution 
of 


Let j be an index that takes values from 1 to 3. Then we have 


T’njtpj + 'Dn4'lp4 + + T^nQlp& = g4>' 5n4—Ipj + g4>' 5i—1p4 ■ 


With the useful relations 


)|, S. 103f.] 

Xj A 1 2 

— +4 = -Cl+j + Co+j 


Ipj = -y+jfi =x{fiYi i+i m + f2YI 1-1 „i + fsYI I m) 
= -flCl+t + f2C0(p\ , 


the right-hand side of (A17) can be written as: 


g(p'Sn4 {-flCl+t + f2Co+t) + g(p'5i {-Cl+] + Co(^f) /4 ■ 


d^ 2 d Injln + 1) 

r dr 


^nk '^nkilP^ 


fk{r) = 0 , 


with V{(p) = 


! m1{(p) 0 0 —cigcp' 0 

0 Wa((/)) 0 Co#' 0 

0 0 ml[(P) 0 0 

-Cl#' Co#' 0 ml{(p,T) 0 

0 0 0 0 -2 

Vo 0 0 0 0 


0 
0 
0 
0 

ml{(p,T) 0 
m 


(A21) 


(A22) 
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These considerations motivate the following ansatz for the Green’s function G{x,x'): 


Gkm{x,x') = ^grs{r,r')ifl{x)ip‘^^{x') , 


(A23) 


where the radial Green’s functions - their /-dependence is not displayed explicitly - satisfy 


92 2 5 /„(/„ +1) 


^ dr'^ r dr 

In (^) one needs to know the following expressions: 

+ °o n, , w 0^2 

47r 1 4 


^nk 


gkm{r, r') = — -5(r - r')(5„ 


tr 




dif) 


G 


— ^ I(511 + 922 + 533) + 2 g 44 


1=0 


+3{Xt 4> — ET)g^5 - 1 - 


tr [c CVG] = ^ XI Cm Vr-gabV^mV’n* ^ ^ ^ ^9]iE^ 


a .J>* _ ^ / V7 „ .G,A* _1_ 

l,m 

3^^^2^ + l^ /, 2 \ ( , 2 


=ii:- 


27r 


— Cl ( 314 + -314 I + Co I 324 + “324 


(A24) 


(A25) 


(A26) 


Obviously the components 333 and 366 satisfy the same differential equation (see ( A24 )). Hence, they can be combined. 
Due to the different sign of their factors of degeneracy they cancel exactly, except for the s-wave. As mentioned above 
at / = 0 the amplitudes f 2 (r) and fair) vanish and therefore 333 = 0 in this particular partial wave. Hence, the 
(6 X 6)- reduces to a (5 x 5)-system, where the fifth channel only contributes at / = 0. The final result is given in 
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APPENDIX: TABLES 


UlH 

'T'h.t 

^ c 


^ c 

30 

40 

67.572 

74.811 

68.238 

75.464 

67.607 

74.916 

68.070 

75.433 


Table 1: The values of the roll-over temperature Tq and the critical temperature Tc obtained from the high 
temperature potential Vht and the approximation of the one-loop resummed effective potential Vr. The 
top-quark mass is chosen as rrit = 170GeV. The Higgs boson mass mn and the temperatures are given in 
GeV. 


APPENDIX: FIGURE CAPTIONS 

Figure 1: The expansion of the one-loop correction 5'^^ to the effective action. The lines represent the propagators 
and the dots the vertex factors given by Wi((/), T) — Ui{0, T). 

Figure 2: Gomparison of the logarithm of the nucleation rate 7 (see (^)), i. e. including the one-loop corrections 
and prefactors, versus the temperature T. 7 has dimension GeV"* and is calculated in units of {gv{T))^. The squares 
represent the results we obtain using the profile (jjci, which is a solution of ([^. They are connected by straight dashed 
lines. The dots are the results we obtain using the the self-consistent solution of the full one-loop field equation (|^). 
The solid line is the function F{T) as defined in (^5|), with a = 3.48 • 10“^ and b = 1.350. The parameter set used 
here is mn = 30GeV and mt=170GeV. 

Figure 3: The same as in figure 2 for the parameter set rriH = 40GeV and TOt=170GeV. Here the solid line is the 
function F{T) as defined in (^), with a = 2.38 ■ 10“^ and b = 1.362. 

Figure 4: Gomparison of the logarithm of the nucleation rate after mapping of the temperature intervals [T[^,rJ^] 
and [Tq*^,T^*^] to [0,1]. For the parameter set niH = 30GeV (to// = 40GeV) and toj = 170GeV the squares (triangles) 
are the results obtained when using the profile (j>c\ and the dots (diamonds) are due the self-consistent profile. 

Figure 5: The logarithm of the fraction pa{t) of space that has been converted to the asymmetric phase in dependence 
of the time t. The three lines are due to three different values of the bubble wall velocity v. The solid line is obtained 
for V = c, the short-dashed line for u = c/10 and the long-dashed one for v = c/1000. The parameter set used here is 
rriH = 30GeV and rrit = 170GeV. 

Figure 6 : The total number of bubbles per unit volume N{t) in dependence of the time t. The three lines are due 
to three different values of the bubble wall velocity v. The description is the same as in figure 5. 

Figure 7: The average bubble radius R{t) in dependence of the time t. The three lines are due to three different 
values of the bubble wall velocity v. The description is the same as in figure 5. 
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